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ABSTRACT 

We propose a novel method to constrain the spatial extent of dust around galax¬ 
ies through the measurement of dust temperature. Our method combines the dust 
emission of galaxies from far-infrared (FIR) image stacking analysis and the quasar 
reddening due to the dust absorption around galaxies. As a specific application of our 
method, we use the stacked FIR emission profiles of SDSS photometric galaxies over 
the IRAS 100/im map, and the recent measurement of the SDSS galaxy-quasar cross¬ 
correlation. If we adopt a single-temperature dust model, the resulting temperature is 
around 18K, which is consistent with a typical dust temperature for a central part of 
galaxies. If we assume an additional dust component with much lower temperature, 
the current data imply the temperature of the galactic dust needs to be higher, 20K to 
30K. Since the model of the density and temperature distribution of dust adopted in 
the current paper is very simple, we cannot draw any strong conclusion at this point. 
Nevertheless our novel method with the elaborated theoretical model and multi-band 
measurement of dust will offer an interesting constraint on the statistical nature of 
galactic dust. 
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1 INTRODUCTION 


mated by a single power-law: 


Dust plays important roles in cosmic star formation and 
evolution of the galaxies. The basic ingredients of dust 
grains are metals produced through past stellar activity, and 
thus the main reservoir of dust is conventionally thought 
to be main ly confi ned in interstellar space within galax¬ 
ies. IZwickv I (119621 1 . however, suggested the existence of 
dust filling the intracluster space within the Coma cluster. 
This motivated the investigation of the abundance and spa¬ 
tial distribution of dust in different environments, including 
the color-excess of backgro u nd ob j ects due to dust optical- 
UV reddening dZaritslY 1994: Chel ouche et al. I 120071 : 

E 


iMcGee fe Balogh I l20ld : iMuller et al. 112 0081), a nd th e FIR 


Stickel et al. 20021: Kaneda et al. 

2001)1: Kitavama et al. 

2009). and from stacking analysis ( 

Montier & Giard 20051: 

Gutierrez & LoDez-Corredoira 20141). 


Recently, Menard et al. (2010a: hereafter MSFR) inves¬ 
tigated the distribution of dust around galaxies by measur¬ 
ing the angular correlation between SDSS galaxy distribu¬ 
tion and distant quasar colors. They found that the mean 
g-i reddening profile around SDSS galaxies is well approxi- 


/ a \ — 0 . 86 ± 0.19 

(S 9 - i )(0) = (1.5±O.4)xlO- 3 (^—J , (1) 

where 6 is the angular separation between foreground galax¬ 
ies and background quasars. Furthermore they discovered 
that the above power-law extends even for 6 > 10\ The an¬ 
gular scale corresponds to several Mpc at the mean redshift 
( z ) = 0.36 of their SDSS galaxy sample. This is far beyond 
the typical scale of galactic disks, and even larger than the 
virial radius of typical galaxy clusters. 

MSFR appear to interpret their result as an evidence 
for an extended dust surrounding an individual galaxy be¬ 
yond a few Mpc, which we refer to as the circum-galactic 
dust model (CGD model). Their interpretation, however, is 
rather subtle. The mean reddening profile from their mea¬ 
surement oc 8~ 0 ' 8 is close to that of the angular correlation 
function of galaxies. Thus the detected dust reddening may 
be equally explained by the summation of the dust compo¬ 
nent associated with the central part of galaxies according 
to the spatial clustering of those galaxies, which will be re¬ 
ferred to as the inter-stellar dust model (ISD model). 

In practice, it is difficult to distinguish between the 
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CGD and ISD models on the basis of the statistical cor¬ 
relation analysis alone as performed by MSFR.. Therefore 
a complementary and independent method to constrain the 
nature of the dust is needed. This is exactly what we attempt 
to propose in this paper. 

For that purpose, we measure the dust far-infrared 
(FIR) emission of the SDSS galaxies by image stacking 
analy sis. Similar analysis on the SFD Galactic extinction 
map dSchlegel. Finkbeiner fe Davislll998l. SFD1 has det ected 

the FIR emission of SDSS galaxies ( Kashiwagi et al]|2013l . 

KYS13). We return to the 100pm intensity map by SFD, 
instead of their extinction map, and perform the stacking 
analysis of the same galaxy sample used by MSFR. If the 
detected FIR emission originates from the same dust com¬ 
ponent as the MSFR reddening measurement, the emission 
to absorption ratio puts a constraint on dust temperature, 
which would in turn offer complementary information to 
distinguish between the CGD and ISD models mentioned 
above. 

The present paper is organized as follows. The data used 
in the current analysis are described in Section 2. In Section 
3, we perform the stacking analysis of the MSFR galaxy sam¬ 
ple on IRAS/SFD 100pm map. We show the constraint on 
the dust temperature from the detected FIR emission com¬ 
bined with the MSFR reddening measurement. We present 
summary and conclusions of the paper, and discuss future 
outlook in Section 4. Throughout the analysis, we assume 
the standard ACDM cosmology with D m = 0.3, Da = 0.7, 
and h = 0.7. 


2 DATA 


We select our galaxy sample from the SDSS DR7 pho¬ 
tometric galaxies with 5 passbands, u, g, r, i, and 
z , in northern galactic cap, which cov ers ^7600 deg 2 


For details of the photometric data, see Stoughton et al 

B - —’ —’ 


t hotom e tnc d a ta, see | bt 
dl998l. 200611 : Fukueita et alJ (1991); 
Hogg__et_al] ( 200l|) ; |lvezic_et_aL ~|2004j) ; Snuth_etaL (2002) 


Tucker et al. ( 2006l lT lPadmanablian et al] ( 20081 1 : lPier et al] 

( 2003h . We conservatively masked ~ 5% of the total area 
following the SDSS mask definition. We also removed the 
objects with bad photometry or fast-moving flag according 
to the photometry flags, which are su spicious to be corre¬ 
lated with the Galactic foreground. See lYahata et all (l2007l l 
for more details of our data selection. 

For the current analysis, we impose the same i-band 
magnitude cut, 17 < rrii < 21, as the MSFR sample for 
a direct comparison with their results, where the magni¬ 
tudes of the galax ies are correct for Galactic extincti on us¬ 
ing the SFD map (ISchlegel. Finkbeiner fe DavisHl998h . Our 
final sample collects 2.88 x 10 7 galaxies. 

For far-infrared data, we use the all-sky diffuse 100pm 
map provided by SFD; they have carefully processed the 
original IRAS/ISSA 100pm sky map, removing the scan 
pattern of IRAS, correcting calibration errors based on 
COBE/DIRBE data, and subtracting zodiacal dust emis¬ 
sion and bright point sources with / 60 /jm > 0.6Jy. 

Hereafter, we adopt a Gaussian with a = 3'T for the 
point spread function (PSF) of SFD/IRAS map, as mea¬ 
sured by similar stacking analysis by KYS13. 


3 IMAGE STACKING ANALYSIS OF FIR 

EMISSION FROM SDSS GALAXIES 

3.1 Stacked radial profiles 

Following the procedures of KYS13, we stack the SFD/IRAS 
100pm map over 120' x 120' squares centered on each SDSS 
galaxy. Each image is randomly rotated around the center. 
The resulting stacked image shows clear circular signature 
of dust emission associated with those galaxies (KYS13). 

The radial profile of the raw stacked image is shown 
in Figure [lj>. The quoted error bars reflect the rms in each 
radial bin (A# = l'.O). The radial profile is reasonably well 
fitted by Gaussian corresponding to the PSF around the 
central region, but exhibits an extended tail beyond the PSF 
width, a = 3 , .l, which corresponds to roughly IMpc for the 
mean redshift (z) ~ 0.36 of the SDSS galaxies. 

At sufficiently large 9, the stacked flux should be dom¬ 
inated by the Galactic foreground, which is uncorrelated 
with the SDSS galaxies and expected to be constant. The 
stacked flux, however, increases beyond 9 > 30'. While we 
do not completely understand the behavior (see also discus¬ 
sion in KYS13), it may be partly due to the fact that the 
SDSS survey region is designed to be located at the low- 
extinction region, therefore towards high galactic latitudes. 
Thus the outskirt of the SDSS region is surrounded by low 
galactic latitudes with relatively higher values of the 100pm 
intensity, and the stacked flux centered at the SDSS region 
tends to be systematically larger at larger 9. Nevertheless 
the profile for 9 < 20' matches nicely that expected from 
angular correlation functions of SDSS galaxies (Okabe et al. 
in preparation). This is why we adopt the profile modeling 
discussed below. 

We adopt the following radial density profile of dust: 

I tot {9) = I s (9) + I c {9) + C, (2) 

where I s and 7 C represent the contributions from the central 
single galaxy (single term) and from the clustered neighbor 
galaxies (clustering term), respectively, and C is the back¬ 
ground level of the foreground Galactic dust emission Q. We 
assume that the Galactic foreground, C, should be uncor¬ 
related with the SDSS galaxies, and thus is assumed to be 
constant at 9 < 20'. 

Since the PSF of SFD/IRAS map is well approximated 
by Gaussian, I s (9) is written as 

W) = Iso exp - ( 3 ) 

where a = 3'.1 is the Gaussian width of PSF. 

The clustering term 7 C is written in terms of 7 S and an¬ 
gular two-point correlation function (2PCF) of galaxy, w(9), 
as 

m = [ d m' d -^l 

J dm 

x / dipl s (0 — y\m)w((p\m), (4) 


1 These definitions of /tot, Is , and I c are equivalent to £g 0t , £|, 
and Eg used in KYS13, respectively, except that E g denotes the 
SFD map extinction in units of [mag], whereas I in this paper 
denotes the intensity in units of [MJr/sr]. 
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where dN(m')/dm' is the differential number count of galax¬ 
ies (whether or not detected by SDSS) as a function of m!. 
We assume that the single term is written as a function of nm 
alone, therefore the dependence on other physical quantities 
is neglected. We approximate the angular 2P CF is described 
as a single power-law in t his angular scale (I Connolly et al.l 
120021 : IScranton et al.|[ 2002 f ): 

w(9;m) = A(m) (J^ , (5) 

where the amplitude A is a function of rrn, but the index 7 
is assumed to be a constant and independent of rrn. In this 
case, equation © reduces to 

«»> = '.»“■>(-•&) < 6 > 

where l F\{a\b\c) denotes the confluent hypergeometric 
function, and 



x J dm 1 I s o(m')A(m') dN J^ - . (7) 

We fit the radial profile of the stacked image using equa¬ 
tions ©, ©, and ©. In doing so, we do not use equation 
0 , but treat I c 0 simply as one of the fitting parameters em¬ 
pirically determined from the observed profile. Consistency 
of the resulting 7 C 0 with equations © and 0 independently 
measured for SDSS galaxies is an interesting topic (KYS13), 
which will be discussed in detail elsewhere (Okabe et al. in 
preparation). We estimate the statistical errors using the 
jackknife resampling method by dividing the entire SDSS 
sky area into 400 patches of equal area. 

The detected emission profile at small 9 is affected due 
to the IRAS PSF, and should not be directly compared with 
the MSFR measurement. Therefore, we use the clustering 
term, which is relevant for 9 a, for the dust temperature 
constraint in the following Section. In fact, the PSF effect 
on the clustering term vanishes at large 9 and equation © 
reduces to the power-law as 

m = r(1-7/2) (7^) ' (8) 


Figure 1. Radial profile of the FIR stacked image of SDSS galax¬ 
ies. The symbols indicate the radial average of the stacked image 
and the error bars show rms in each radial bin. (a); Radial pro¬ 
file of the raw stacked image for 0 < 60' (before subtracting the 
offset level due to the Galactic dust). Shaded region indicates 
the data with 6 > 20 / that are not used in the current analysis. 
( 6 ); Radial profile of the stacked image after subtracting the off¬ 
set level of C = 1.39927[MJy sr —1 ], which is computed assuming 
7 = 0.86 and shown as the dashed line in panel a. The lines indi¬ 
cate the best-fits for I s (black dotted), I c (red solid), and /tot — C 
(blue dot-dashed); see equation j2j- Just for reference, we plot 
the power-laws of 7 = 0.65, 0.86 and 1.05 in dashed lines, which 
covers the ranges of the MSFR result in equation q. 


Since we (implicitly) assume here that the mean red¬ 
dening profile of MSFR, equation 0 , is explained in the 
clustered dust model, the value of 7 in equation © should 
match the MSFR result. In order to confirm the validity of 
the assumption, we first choose 7 s o, 7 c o, C , and 7 as free 
parameters, and fit to the observed profile imposing I s 0 ^ 0 
and Ico ^ 0. The resulting best-fit value, 7 = 1.07 ± 0.16, 
is consistent with that of MSFR, 7 = 0.86 ± 0.19 (the other 
best-fit values include 7 s o = 0[MJy sr -1 ], 7 c o = (3.5 ±0.4) x 
10 -3 [MJy sr -1 ], and C = 1.399 ± 0.035[MJy sr -1 ]). Indeed 
as Figure [T] illustrates, the difference among the predicted 
profiles for 0.65 < 7 < 1.05 is very small for the angular 
scales of our interest 9 a. The departure from the power- 
law for I c < 6 x 10 -4 [MJy sr -1 ] is not a problem because it 
simply reflects the sensitivity to the subtracted offset C. 

Thus we fix 7 = 0.86 in what follows, and obtain the 
best-fit parameters as 7 a o = (6.1 ± 4.0) x 10 -4 [MJy sr -1 ], 
Ico = (3.1 ± 0.7) x 10 -3 [MJy sr -1 ], and C = 1.399 ± 
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0.035[MJy sr -1 ]. The best-fit profile for each component is 
shown in Figure [T] 

The stacked FIR emission profile corresponding to the 
clustering term for 7 = 0.86 is finally given as 


(7l00 M m)(#) 


(7.0 ± 1.6) x 10 ~ 3 / 6 \~ 0 86 
MJy sr _1 Vl'.O / 


(9) 


at large 8, which plays a major role in our method proposed 
in Section EU 

We note that while the statistical error of C is much 
larger than the best-fit values of 7 s o and 7 c o themselves, it 
does not affect the detection significance of the dust emission 
from SDSS galaxies. In fact, the variance of C simply comes 
from that of the Galactic dust over the SDSS survey area; 
the majority of the 400 jackknife subsamples indicates simi¬ 
lar signatures of the dust emission, except for the difference 
of C. 

It is interesting to note that the observed stacked profile 
and the prediction from the summation of individual SDSS 
galaxies indeed agree very well, as mentioned in Section 1. 
Thus we will proceed further to an independent and com¬ 
plementary analysis in order to constrain the spatial extent 
of dust in the rest of this section. 


3.2 A simple model prediction of the dust 
emission 


While the dust extinction is determined mainly by its col¬ 
umn density, the dust emission depends sensitively on its 
temperature as well. Therefore, if the measured extinction 
and emission comes from the same dust distribution, their 
ratio serves as a sensitive measure of the dust temperature. 
In this subsection, we will explicitly show theoretical ex¬ 
pressions for the reddening and emission of dust in a simple 
model of dust density distribution. Since we are interested 
in the scales beyond the galactic disk scale, we consider the 
clustering term alone. 

The angular profiles of dust extinction and emission 
around a galaxy are calculated by integrating the dust sur¬ 
face density Ed (r p ,z) of nearby galaxies at 2 separated by 
the projected distance r p = (Ia(z)9 from the central galaxy, 
where cIa(z) is the angular diameter distance at z. For sim¬ 
plicity, we assume that the 2 -dimensional projected dust sur¬ 
face density responsible for the clustering term is given by 
a single power-law: 


E d (r p ,z) = E d 0 (z) \ (10) 

\ r p, 0 / 

Throughout the current model, we set the power-law index 
7 as that of the galaxy angular correlation function, equa¬ 
tion ©, specifically 7 = 0.86 in what follows. Although 
we neglect the redshift evolution of the correlation length 
r P: o, it is effectively absorbed in E d o(z) as long as 7 is time- 
independent as assumed here. 

Under the above assumptions, the angular extinction 
profile of dust at redshift z is written as 


E g -i(e,z) = 


x 


2.5 
In 10 
2.5 
In 10 


T 


x 3 \ 

IT z) 




lTz 


E d(d,A(z)8, z), 


-rle 


Xi 

lTz 


^ext 



(ii) 


where X g and A; are the rest-frame wavelengths of SDSS 
g and i-bands, respectively, and K ex t(A) is the extinction 
cross-section per unit dust mass at a wavelength of A. The 
average angular extinction profile around SDSS galaxies is 
then given by 


(E g -i)(9) = 


x 


2.5 
In 10 



[/" 

.lyo 

K-ext 


x E d o(z) 


dN , ' 

dz 

dz 

' Xg \ 
l + z) 
dA{z)d 
r P ,o 


1 


Kext (irl) 

7 dN , 

——dz, 


( 12 ) 


where dN/dz is the redshift distribution of SDSS 
galaxies. Following M SFR, we adopt an approximation 
(IDodelson et al][2002l k 


dN 

dz 


oc z 2 e -(*/°- 187)1 ' 26 . 


(13) 


Thus the number-weighted mean redshift of the sample is 
given by 


<*> 


f z(dN/dz)dz 
f ( dN/dz)dz 


(14) 


One can similarly compute the angular FIR emission 
profile around SDSS galaxies. Since the dust emission at A = 
100 /rm is well approximated by the blackbody spectrum, the 
corresponding surface brightness at redshift z is given as 


7l00 P m($, z 1 T d ) 


1 

0TW 


B v 


(100 /xm ^ 

ItttW 


x «.b. (x/T) Sdp/iMM, (is) 


where K a b s is the absorption cross section per unit dust mass, 
7?„ is the blackbody spectrum per unit frequency, T d is the 
dust temperature, which we assume to be independent of z, 
and the same for all SDSS galaxies, and 1/(1 T z ) 4 comes 
from the cosmological dimming effect. 

The average angular emission profile of SDSS galaxies, 
which corresponds to I c {9) observed by the stacking analy¬ 
sis, is given as 


(dl00pm)(8, Td) 

dN , ' 

L 

f lOO/rmA 

Kabs [lT7) 


f 

Sdo(^) 


(1 + z ) 4 Bl 

dA(z)8 

r P , 0 


f 100/rm 

V 1 + 2 ’ 

~ 7 dN , 

—— dz. 

dz 


T d 


(16) 


Because we adopt the power-law dust profile, equation 
( 11011 . the ratio of equation (1161) to ( 1121 ) is independent of 8, 
and written in terms of K e xt, Kabs, and T d alone. 

The observed profile of the emission to reddening ra¬ 
tio is shown in Figure [5] The filled circles are plotted using 
the residual of the emission profile, from which the best- 
fit single term and the offset level assuming 7 = 0.86 are 
subtracted. The red solid curve shows the ratio of the best- 
fit clustering term, 7 C (0), with 7 = 0.86 to equation ([I|). 
and the shaded region indicates its uncertainty due to the 
statistical error of 7 c o and the amplitude of equation ©. 
The uncertainty of the power-law index in equation © is 
not considered here. At small 9, the emission profile is sup¬ 
pressed due to the SFD/IRAS PSF effect, whereas the ratio 
converges to a constant at large scale. The emission to red¬ 
dening ratio at large 8 limits is given by equations 0 and 
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9 [arcmin] 


dust temperature: Td [K] 


Figure 2. Radial profile of 100 fim emission to reddening ratio. 
Filled circles (black) indicate the observed 100/zm emission profile 
of the stacking analysis divided by the MSFR reddening profile, 
where the best-fit single term I s ( 0 ) and the offset level C assuming 
the 2PCF index as 7 = 0.86 are subtracted. Red solid curve is the 
best-fit profile of the clustering term / c ( 0 ) with 7 = 0 . 86 , divided 
by Eg—i(6) of MSFR. Shaded region indicates the uncertainty 
of the emission to reddening ratio where the statistical error of 
I c 0 and the MSFR measurement is taken into account. Crosses 
(filled triangles) and dashed (dot-dashed) curve indicate the same 
as filled circles and solid curve, but for 7 = 0.65 (1.05). 


© as (Iiooiim}/(Eg-i) = 4 -7± 1-6 [MJy sr _ 1 mag _1 ], which 
corresponds to the shaded regions in Figure [3] below. 

We also consider the extent to which this result is sen¬ 
sitive to the choice of the power-law index 7 , which is fixed 
as 0.86 in the analysis above. We repeat both the fitting 
to the observed profile and the theoretical calculation of 
equation m and (Hhd, varying the value of 7 from 0.65 
to 1.05. Figure [2] shows the observed emission to reddening 
ratio for 7 = 0.65, 0.86 and 1.05. The average ratio changes 
approximately ~ 20 per cent (and its fractional uncertainty 
is similar to that for the case of 7 = 0 . 86 , although it is not 
shown in Figure©. We also make sure that the theoretical 
value from equation m and m changes by 10 per cent 
according to the corresponding change of 7 . Consequently, 
we find that the uncertainty of dust temperature due to the 
choice of 7 is merely ~ IK. This can be neglected compar¬ 
ing with the possible larger systematics due to other many 
simplifying assumptions. 


3.3 Constraints on dust temperature 

The solid and dashed lines in Figure © indicate the ex¬ 
pected emission to extinction ratio as a function of Td- We 
adopt the values of K ex t and K a bs, from the dust model by 


Figure 3. Constraints on the FIR emission to extinction ratio 
from MSFR and the stacking results (shaded region). Solid and 
dashed lines indicate the prediction for Mi lk y Wa y (Ry = 3.1) 
and SMC dust model jWcingartner fc Draine r200 lli , respectively. 
The power-law index of galaxy 2PCF is fixed as 7 = 0.86. If 
7 = 0.65 and 1.05 is assumed, the allowed region is shifted by 
—20 and 20 per cent, respectively. 


IWeingartner fc Draine I d 2001 1 ) 0 for Milky Way (Rv = 3.1) 
and SMC dust, for solid and dashed lines, respectively. 

Here the redshift dependence of Edo(z) is neglected and 
assumed to be constant just for simplicity. Indeed we made 
sure that the ^-dependence of Edo(-z) does not significantly 
change the result; if we assume Edo(a) oc (1 + z) p for in¬ 
stance, the model prediction of ( 1 10011m}/(E g -i) changes by 
q=15 per cent for p = ± 1 , and the dust temperature con¬ 
straint changes by T0.2K . The r ecent measurement of dust 
mass function by iDunne et al.l (1201 ll ) found that the cos¬ 
mic dust mass density in sub-mm galaxies rapidly increases 
with redshift up to z ~ 0.5. Thus the constraint on the dust 
temperature below may be slightly underestimated. 

Figure ©indicates that the dust model predictions and 
the observed region are consistent if Td = 18.2j)o'gK for MW 
dust, and Td = 17.0j7gK for SMC, thus the obtained con¬ 
straints are almost insensitive to the choice of dust model. 

Given several approximations adopted in our simple 
model, the quoted statistical errors may underestimate the 
real uncertainty of the dust temperature. Nevertheless it is 
encouraging that the derived dust temperature is in good 
agre ement with that of the typical cold component of the 
ISP (ISc hlegel. Finkbein er fe Davislll998l : IDunne et al.ll20fll : 
Idemens et al. 2013h . While this result is reasonably con¬ 
sistent with the ISD model, it is premature to conclude that 
the observed dust profile is explained by the sum of the dust 
in the central parts of individual galaxies. Our dust model 


2 Data is taken from Web-site of B. T. Draine, 
http: //www.astro, princeton.edu/~draine/dust/dustmix.html 
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above is based on a single-temperature component, and the 
extended dust component around a galaxy may have a sub¬ 
stantially lower temperature. In this case, the emission flux 
is much smaller while the reddening amplitude remains al¬ 
most the same. 

Therefore we consider a two-component dust model be¬ 
low. In order to keep the same surface density profile of dust, 
we assume that the two components, corresponding to ISD 
and CGD, share the identical spatial distribution, but they 
are allowed to have different temperatures Tisd an d Tcgd- 
While this is still a very simple model, we would like to 
proceed with it because the low-angular resolution of IRAS 
images make it difficult to distingush the density profile less 
than a = 3'.1. In any case the purpose of the current paper 
is to propose a new method to constrain the dust density 
and temperature, which will be improved significantly later 
theoretically and observationally. 

In the above spirit, we replace equation (EH by 


Sd(r p , z) = [£dO,ISD(z) + SdO.CGD^)] 


r P , o 


(17) 


Then the observed FIR emission to extinction ratio becomes 

(flOO/jm) 


( Eg-i) \ oba 

,/"lSD (7l00 P m) (Tisd) + (1 — /lSD)(7l00Mm)(TcGD) 
{Eg-i) 


(18) 


at large 6, where (E g -i) and (/iooijm)(T) in the right-hand- 
side are given by equation m and m, respectively. We 
further assume that the fraction of the ISD mass to the total 


/isd = 


SdO.ISD 


SdO.ISD + SdO.CGD 


(19) 


is independent of redshift. Thus the observed value of 
(Iiooiim)/{Eg-i} provides a constraint on /isd for given val¬ 
ues of Tisd and Tcgd- The preceding analysis corresponds 
to /isd = 1. 

The temperature of the CGD is fairly uncertain be¬ 
cause of the unknown heating mechanism of the CGD. If 
the heating source of the CGD is dominated by the cos¬ 
mic UV background, which is lower than the interstel¬ 
lar radiation field in the solar neighborhood by two or¬ 
ders o f magnitude dMadau fe Pozzettill200CI : iGardner et al.l 
l2000l: IXu et ali l200i) ~ we obtain Tcgd = 10K following 
jDrain^^Le^JlflSjMDrain^ 1201 il l. On the other hand, 
lYamada fe Kitavama 1 1 2005ll assume the collisional heating 
mechanism by hot plasma and the efficient injection of dust 
grains outside the galactic disk. In this case, they suggested 
a possibility that the dust temperature reaches even ~ 30K. 

The abundance of such possible high-temparature 
CGD, however, is s everly constrained by the observed data. 
(tClark et al. Il2015l l fitted the SED of dust-selected galaxies, 
and found that the majority of them have cool and warm 
dust components with T ~ (10-20)K and > 30K. The mass 
of the cool component is typically 100 times larger than 
that of the warm component. While the two components 
may correspond to the two-temperature phases in ISD, the 
estimated mass ratio can be also interpreted to put a severe 
constraint on the presence o f the hot CGD with T ~ 30K. 
Furthermore, iDraine et al J ( 2014 1 reported that the dust 
temperature near the edge of M31 disk is 15K. Thus the 



Figure 4. Constraints on the ISD temperature, Tisdi and the 
mass fraction of the ISD to the sum of the ISD and CGD. We 
assume that the CGD temperature as Tcgd = 10 K. Red curves 
and blue shaded region indicate the co nstraint adopting Milky 
Way (Ry = 3.1) and SMC dust model dWeingartner fe DraincJ 
|2001 : 1. respectively. 


dust temperature in the outskirt is naturally expected to be 
much lower. 

For those reasons, we adopt Tcgd = 10K in what fol¬ 
lows. We also confirm that the result below does not change 
as long as Tcgd is lower than 10K. 

Figure[4] shows the constraint on Tisd and /isd from the 
observed value of (IiooMm)/(-E 9 -i). Due to the strong degen¬ 
eracy between the two parameters, the constraint allows a 
wide range of the ISD mass fraction, even as small as /isd ~ 
10~ 2 if Tisd = 30K. Thus the measurement of the mean dust 
temperature of the central parts of galaxies Tisd is crucial 
in distinguishing the origin of the spatial extension of dust. 
Indeed, the average temperature of the ISD varies depend¬ 
ing on the properties of the g alaxies, and can be as high 
as ~ 40K (ISkibba et al. 1120111 ). In this respect, the current 
data do not exclude a possibility that a substantial amount 
of the CGD exists, a s suggested b y MSFR . and the sub¬ 
sequent studies (e.g., Fukugita|_ 2011 1 _ Menard fe Fukueita I 


120121 : IPeek. Menard, fe Corrales 1 12014 1 . Nevertheless fur¬ 
ther improvements in model predictions and the observa¬ 
tions in future would put more stringent constraints on the 
spatial extent of dust through the measurement of dust tem¬ 
perature as we proposed in this paper. 


4 SUMMARY & CONCLUSIONS 

The spatial distribution of dust is of fundamental impor¬ 
tance in understanding the star formation and metal circu¬ 
lation history in the universe. It is also crucial in correcting 
for the magnitude of distant objects due to the resulting 
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reddening/e xtinction dAguirre 11 19991 : IMenard et al.l l2010bl : 
iFang et al. Il201lh . 

In a previous paper (iKashiwagi et all 120131 / , we have 
detected the FIR dust emission from SDSS galaxies via their 
image stacking analysis. We found that the amount of dust 
emission is largely responsible for the observed anomaly in 
the surface den sity of SDSS galax i es as a function of th e 
SFD extinction (jYahata et alJl2007l : iKashiwagi et al]l201fH b 
Our previous analysis implicitly assumed that the dust of 
each galaxy is locally confined in the galactic disk scale, and 
that the observed FIR emission within the large PSF width 
(FWHM= 6'.1) is simply given as a sum of contributions of 
individual galaxies (corresponding to the ISD model in the 
present paper). In contrast, the dust around a galaxy may 
be indeed spatiall y exten ded up to ~ IMpc (CGD model), 
as claimed b y IM enard et al.l (12 01081 1 and more recently by 
IPeek. Menard, fc Corrales I (120141 / through the correlation 
of background object colors against the separation length of 
foreground galaxies. 

In order to distinguish between the ISD and CGD mod¬ 
els, we propose a new method that constrains the tempera¬ 
ture of dust by combining the absorption (detected through 
reddening of quasars) and emission (detected through the 
stacking of galaxies) features. Assuming that the nature of 
galactic dust is described by those of MW and SMC, we 
find that the observed dust is reasonably explained in terms 
of ISD model if the dust temperature of the central parts 
of individual galaxies, Tisd, is approximately 20K. The es¬ 
timated temperature is consistent with that of the galactic 
dust in the central region, but may be higher than that pre¬ 
dicted for CGD, if it is heated by UV background alone. On 
the other hand, the substantial amount of dust may reside 
far outside the galactic disks if Tisd is much higher than 
20K. 

Given several simplification and approximations that 
we adopted in the present simple model analysis, the asso¬ 
ciated error-bars of the derived dust temperature is fairly 
uncertain. Nevertheless we would like to emphasize that the 
main purpose of the present paper is to propose a new ob¬ 
servational method to diagnose the nature of galactic dust. 
Therefore we do not discuss the interpretation of the present 
preliminary result. 

Our proposed method should be, and indeed can be, 
improved in many ways; the two components of dust may 
have different spatial density profiles in addition to the dif¬ 
ferent temperatures. The redshift evolution of the tempera¬ 
ture and amount of dust may be included in the theoretical 
models. In addition, the line-of-sights of our emission and 
reddening measurements may not be exactly the same. Since 
quasars behind the heavily extincted line-of-sights may not 
be identified in the SDSS photometric catalogue, the red¬ 
dening measurement may systematically underestimate the 
real optical depths while the emission is largely free from 
the bias. Those improvements of the theoretical models and 
the effect of the possible selection bias need to be investi¬ 
gated, for instance, with mock simulations, which is beyond 
the scope of the present paper. 

The observational data and analysis can be also im¬ 
proved in future. The dust temperature would vary depend¬ 
ing on the different properties of galaxies, and the amount 
of dust emission should depend on the morphology of galax¬ 
ies. The image stacking analysis with better angular resolu¬ 


tions and in multi-wavelengths would significantly improve 
the observational data. Indeed current result is significantly 
limited by the poor angular resolution of IRAS. In those 
respects, the higher-angula r-resolution and multi -band far- 
infrared data by AKARI dMurakami et al] 120071 / are very 
promising. We plan to present elsewhere more detailed and 
systematic results using the AKARI data (Okabe et al. in 
preparation). 

We thank Brice Menard, Bruce T. Draine, Masato Shi- 
rasaki, and Tetsu Kitayama for useful discussions. We also 
appreciate several comments by an anonymous referee that 
motivated us to consider the two-temperature dust model 
as well. This work is supported in part from the Grant-in- 
Aid No. 20340041 by the Japan Society for the Promotion 
of Science. Y.S. and T.K. gratefully acknowledge supports 
from the Global Scholars Program of Princeton University 
and from Global Center for Excellence for Physical Science 
Frontier at the University of Tokyo, respectively. 

Data analysis were in part carried out on common 
use data analysis computer system at the Astronomy Data 
Center, ADC, of the National Astronomical Observatory of 
Japan. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In¬ 
stitutions, the National Science Foundation, the U.S. De¬ 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is http://www.sdss.org/ 

The SDSS and SDSS-II are managed by the Astrophys- 
ical Research Consortium for the Participating Institutions. 
The Participating Institutions are the American Museum of 
Natural History, Astrophysical Institute Potsdam, Univer¬ 
sity of Basel, Cambridge University, Case Western Reserve 
University, University of Chicago, Drexel University, Fermi- 
lab, the Institute for Advanced Study, the Japan Partici¬ 
pation Group, Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Particle As¬ 
trophysics and Cosmology, the Korean Scientist Group, the 
Chinese Academy of Sciences (LAMOST), Los Alamos Na¬ 
tional Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), 
New Mexico State University, Ohio State University, Uni¬ 
versity of Pittsburgh, University of Portsmouth, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 


References 

Aguirre, A. N. ApJL, 1999, 512, 19 
Clark, C. J. R., et al. 2015, astro-ph/1502.03843 
Clemens, M. S., et al. MNRAS, 2013, 433, 695 
Chelouche, D., Koester, B. P., & Bowen, D. V. 2007 ApJ, 
671, 97 

Connolly, A. J., et al. 2002, ApJ, 579, 42 
Dodelson, S., et al. 2002, ApJ, 572, 140 
Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 
Draine, B., T. 2011, Physics of the Interstellar and In- 
tergalactic Medium (Princeton, NJ: Princeton University 
Press) 






















8 Kashiwagi & Suto 


Draine, B., T. et al. 2014, ApJ, 780, 172 
Dunne, L., et al. 2011, MNRAS, 417, 1510 
Fang, W., Hui. L., Menard, B., May, M., & Scranton, R. 
2011, PhRvD, 84, 3012 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shi- 
masaku, K. & Schneider, D. P. 1996, AJ, 111, 1748 
Fukugita, M., 2011, astro-ph/1103.4191 
Gardner, J. P., Brown, T. M., & Ferguson, H. C., 2000, 
ApJ, 542, 79 

Gunn, J. E., et al. 1998, AJ, 116, 3040 
Gunn, J. E., et al. 2006, AJ, 131, 2332 
Gutierrez, C. M., & Lopez-Corredoira, M. 2014, A&A, 571, 
66 

Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J. & Gunn, 
J. E. 2001, AJ, 122, 2129 

Ivezic, Z., et al. 2004, Astronomische Nachrichten, 325, 583 
Kaneda, H., Yamagishi, M., Suzuki, T., & Onaka, T. 2009, 
ApJ, 698, 125 

Kashiwagi, T., Yahata, K., & Suto, Y. 2013, PASJ, 65, 43 
(KYS13) 

Kashiwagi, T., Suto, Y., Taruya, A., Kayo, I., Nishimichi, 
T., & Yahata, K. 2015, ApJ, 799, 132 
Kitayama, T., et al. 2009, ApJ, 695,1191 
Madau, P., & Pozzetti, L. 2000, MNRAS, 312, 9 
McGee, S. L., & Balogh, M, L. 2010 MNRAS, 405, 2069 
Menard, B., Scranton, R., Fukugita, M., & Richards, G. 
2010a, MNRAS, 405, 1025 (MSFR) 

Menard, B., Kilbinger. M., & Scranton, R. 2010b, MNRAS, 
406, 1815 

Menard, B. and Fukugita, M. 2012, ApJ, 754, 116 
Montier, L. A., & Giard, M. 2005, A&A, 439, 35 
Muller, S., Wu, S. Y., Hsieh, B. C., Gonzalez, R. A., 
Loinard, L., Yee, H. K. C., & Gladders, M. D. 2008, ApJ, 
680, 975 

Murakami, H., et al. 2007, PASJ, 59, 369 
Padmanabhan, N., et al. 2008, ApJ, 674, 1217 
Peek, J., E., G., Menard, B., & Corrales, L. 2014, 
arXiv:1411.3333 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., 
Kent, S. M., Lupton, R. H., & Ivezic, Z. 2003, AJ, 125, 
1559 

Schlegel, D., Finkbeiner, D., & Davis, M. 1998 AJ, 500, 525 
Scranton, R., et al. 2002, ApJ, 579, 48. 

Skibba, R. A., et al. 2011, ApJ, 738, 89 
Smith, J. A., et al. 2002, AJ, 123, 2121 
Stoughton, C., et al. 2002, AJ, 123, 485 
Stickel, M., Lernke, D., Mattila, K., Haikala, L. K., &Haas, 
M. 1998, A&A, 329, 55 

Stickel, M., Klaas, U., Lernke, D., & Mattila, K. A&A, 
2002, 383, 367 

Tucker, D. L., et al. 2006, Astronomische Nachrichten, 327, 
821 

Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 
Xu, C., K., Donas, J., Arnouts, S., et al. 2005, ApJ, 619, 
11 

Yahata, K., Yonehara, A., Suto, Y., Turner, E.L., Broad- 
hurst, T., & Finkbeiner, D. 2007, PASJ, 59, 205 (Y07) 
Yamada, K., & Kitayama, T. 2005, PASJ, 57, 611 
York, D. G., et al. 2000, AJ, 120, 1579 
Zaritsky, D. 1994, AJ, 108, 1619 

Zwicky, F. 1962, in McVittie G. C., ed., Proc. IAU Symp. 


15, Problems of Extra-Galactic Research. Macmillan, New 
York, p. 347 


